To appear in the Proceedings of Spectroscopies of Novel Superconductors, Cape Cod, September 14-18, 1997 



Andreev Bound States, Surfaces and Subdominant Pairing 
in High Tc Superconductors 
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A characteristic feature of tiie BCS theory of superconductivity is the quantum-mechanical co- 
herence of particle and hole states. Direct observation of particle-hole coherence in unusual su- 
perconducting materials is a strong indication of traditional superconductivity. We use the Fermi 
liquid theory of superconductivity to study the implications of particle-hole coherence on properties 
of d-wave superconductors near surfaces. Typical surface phenomena are the suppression of the 
superconducting order parameter, surface bound states associated with Andreev reflection, anoma- 
lous screening currents, and spontaneous breaking of time-reversal symmetry. We review these 
phenomena and present new results for the effects of surface roughness. 



It is generally accepted that high Tc cuprate supercon- 
ductors are strongly correlated metals whose supercon- 
ducting state is not well described by mean field theory or 
equivalent methods designed for weakly correlated elec- 
trons. On the other hand, theoretical methods for study- 
ing strongly correlated electrons, developed in the con- 
text of high Tc superconductivity, are still far from being 
able to calculate subtle effects such as the influence of sur- 
face roughness on superconducting properties. A theory 
which is not restricted to weakly correlated systems but 
nevertheless can be used to study a broad range of super- 
conducting phenomena including the effects of surfaces is 
the Fermi liquid theory of superconductivity. The most 
useful formulation of this theory is in terms of quasiclas- 
sical transport equations derived in 1968 by Eilenberger 
iQ and Larkin & Ovchinnikov ||^ for superconductors in 
equilibrium, and generalized to nonequilibrium phenom- 
ena a few years later ^,Ql. The Fermi liquid theory of 
superconductivity combines the motion of quasiparticles 
along classical trajectories (external degrees of freedom) 
with the quantum dynamics of internal degrees of freedom 
which are the spin and the particle-hole degrees of free- 
dom. This combination of classical and quantum physics 
is the proper generalization of Landau's semi-classical 
transport equation for normal Fermi liquids to the su- 
perconducting state. Like Landau's theory the leading 
order terms in the expansion parameters of Fermi liquid 
theory are low-energy compared with the Fermi energy 
or long wavelength compared with the atomic scale, e.g. 
1/kf^o, where ^o — hv-f/2'KTc [||J|]- The equations of 
the Fermi liquid theory of superconductivity consist of 
the quasiclassical transport equation for the quasiclassi- 
cal propagator, g(j)f,R; e, i), Eilenberger's normalization 
condition for g, the self-consistency equations for the self- 
energies, a(j)f, R;e,ty, and boundary conditions at sur- 
faces and interfaces |§,0,||Q . An important part of the 
self-energy is the superconducting order parameter A; 
it establishes the coherence between particles and holes. 
We refer to recent reviews |||,g] and the original publica- 
tions 10-0] for the definitions and physical interpretation 



of the quasiclassical propagators and self-energies, and 
the detailed form of the quasiclassical transport equa- 
tions and boundary conditions. 

In quasiclassical theory an excitation approaches the 
surface along a classical (straight) incoming trajectory 
and is refiected into an outgoing trajectory. At a specu- 
lar surface the outgoing trajectory is fixed by the conser- 
vation of parallel momentum (ideal reflection) . Surfaces 
with roughness lead to a statistical distribution of outgo- 
ing trajectories. This classical picture for the kinematics 
of an excitation must be supplemented by the quantum 
equations for the internal degrees of freedom. The in- 
ternal state along a classical trajectory is obtained by 
solving the quasiclassical matrix transport equations on 
this trajectory. The most important quantum effect in 
this context is Andreev reflection [Q, which is caused 
by rotations of the internal state of an excitation from 
particle-type to hole-type (or vice versa). This may lead 
to a velocity reversal (retro-reflection) or to trapping of 
an excitation (Andreev bound states). Particle-hole ro- 
tation and Andreev reflection are controlled by the off- 
diagonal self-energy (order parameter), and carry infor- 
mation about the anisotropy and symmetry of the order 
parameter. A typical trajectory at a (120) surface of a 
dx2_y2 superconductor is shown in Fig.l. For this trajec- 
tory repeated Andreev reflections lead to a bound state 
with excitation energy e = 0, a zero-energy bound state 
(ZBS), which is a robust feature of the spectrum depend- 
ing only on the change in sign of the order parameter 
along the trajectory ||l^ . 

A quasiclassical study of surface effects in unconven- 
tionally paired superfluids (e.g. p-wave pairing in '^He) 
was published by Ambegaokar et al. ||l^], who used de- 
Gennes' method of classical correlation functions |1J] 
to study the suppression and reorientation of the or- 
der parameter. DeGennes' method is a predecessor 
of the full quasiclassical theory applicable to the limit 
A — > 0. The order parameter and the tunneling spec- 
tra at ideal surfaces of a p-wavc superconductor in the 
Balian-Wcrthamer state were calculated by Buchholtz 



and Zwicknagl [|5[. These authors found surface bound 
states and, in particular, a ZBS for trajectories at per- 
pendicular incidence; for the Balian-Werthamer state the 
conditions of the Atiyah-Patodi-Singer theorem are ful- 
filled only at perpendicular incidence. The ZBS shifts for 
all other angles of incidence. 
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FIG. 1. Formation of an Andreev bound state at an ideal 
(120) surface. An excitation moving along the trajectory in 
Fig.la experiences changes in the order parameter A{pf,R), 
where pf is the momentum of the excitation, and R is the 
position. These changes are due to 1) the depletion of the 
order parameter near the surface and 2) the change in mo- 
mentum direction of the excitation when hitting the surface. 
The change in momentum direction leads to a sign change of 
the order parameter as shown in Fig. lb. Fig.lc shows a sketch 
of the order parameter along the trajectory of the excitation. 

The chances for observing zero-energy bound states are 
more favorable for d-wave pairing jl^l. In this case the 
ZBS exists at an ideal surface for a whole range of angles. 
For (110) surfaces and dx2_y2 pairing in bulk this range 
covers all angles of incidence, while the range collapses 



to zero at the (100) and (010) surfaces. These results 
were confirmed and extended by several authors [p7|-fi9[ . 
The ZBS leads to a zero-bias conductance peak which 
was observed in tunneling experiments on cuprate high- 
Tc superconductors by several groups |p0|-p^ . The zero- 
bias conductance peak has been identified convincingly 
as due to tunneling into the ZBS [p5|-p7[. Hence, the 
existence of a ZBS associated with an unconventional d- 
wave order parameter in high-Tc superconductors seems 
well established by experiment as well as theory. 

The order parameter near Tc is determined by the so- 
lution of the linearized gap-equation at Tc. The symme- 
try of this solution defines the dominant pairing channel; 
subdominant pairing channels may mix spontaneously 
in bulk superconductors below a transition temperature 
Tsufc < Tc- However, no such phase transition has been 
observed in high-Tc superconductors. Thus, either there 
is no attractive channel other than the dominant one, 
or the subdominant order parameter is blocked in bulk 
by the dominant order parameter. In the latter case, 
a subdominant superconducting order parameter which 
is not subject to surface pair breaking (s-wave pairing, 
etc.) may be stabilized at surfaces where the dominant 
order parameter is suppressed. This possibility has stim- 
ulated theoretical work by Matsumoto and Shiba ||l^ and 
Buchholtz et al. p9[ on the mixing of subdominant pair- 
ing channels at surfaces. Of particular interest is the 
possibility of time-reversal symmetry breaking by the 
subdominant order parameter, as first discussed in the 
framework of Ginzburg-Landau theory by Sigrist et al. 
p8[ This effect was studied in the full temperature range 
by Matsumoto and Shiba [fsl who used the quasiclassical 
theory to calculate the structure of the order parameter, 
the quasiparticle excitation spectrum at the surface, and 
the spontaneous surface currents. 
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FIG. 2. Calculated phase diagram for the surfaces phases of 
type d^2_y2 +is and d^2_y-2 +idxy at an ideal (110) surface of 
a d^ 
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_y2-superconductor with a bulk transition temperature 
Tc- The bulk transition lines are also shown. 



Fig. 2 shows the phase diagram for the transition to a 
time-reversal breaking sm'face phase at the (110) surface 
of a da;2^u2 -superconductor, as calculated by Fogelstrom 
et al. |29|. The transition temperatures of the surface 



phases are given as a function of the strength of the sub- 
dominant pairing interaction, which we measure by the 
subdominant transition temperature Tc2- Also shown 
are the corresponding bulk transitions. We note that the 
order parameters with s or d.j;y symmetry are not sup- 
pressed by an ideal (110) surface, whereas the bulk stable 
phase, of type dr^2_y2, is suppressed by the surface. As a 
result the sub-dominant order parameter can nucleate in 
a region of strong pair-breaking, and the corresponding 
surface phase spontaneously breaks time-reversal symme- 
try for T <Ts. 



The tunneling density of states at a (110) surface is 
shown in Fig. 3. This figure demonstrates the charac- 
teristic splitting of the ZBS due to time-reversal break- 
ing surface phases, and the effects of surface roughness. 
Note the presence of two pairs of time-reversed Andreev 
bound states for the dx2_y2 +idxy surface phase. Surface 
roughness is described by Ovchinnikov's model, m The 
degree of surface roughness is measured by the parame- 
ter p Q, where the ideal surface corresponds to p = 0. 
The state dx2_y2 + id^y is clearly more sensitive to sur- 
face roughness than the state dx2_y2 -f is. An important 
consequence of time-reversal breaking surface states is a 
spontaneous surface current which flows within a depth 
of a few ^0- The two members of a pair of time reversed 
states, dj.2_y2 + is and dx2_y2 — is or dx^-y2 -\- idxy and 
dx2-y2 — idxy, lead to surface currents of opposite di- 
rection. Typical results for the current density at a flat 
surface, obtained from Fermi-liquid theory of supercon- 
ductivity, are shown in Fig. 4. 




-0.5 0.0 0.5 

e/A„ 



1.5 



0.6 



0.4 



0.2 



0.0 





d+is 
T.,=0.3T 


r 




p=0.0 

p=0.2 

p=0.5 
p=0.8 


//■■- 




•i!5nii!5""^^-rv,i;i 


^/:: 



-10 



-5 




0.06 



0.04 



K 0.02 



0.00 



-0.02 



r =0.37 

c2 L 



-10 




p=0.00 
p=0.02 
p=0.05 
p=0.08 
p=0.20 



FIG. 3. Tunneling density of states at a (110) surface as 
function of energy (e). Numerical results are presented for 
Tc2/Tc = 0.3, T/Tc = 0.1, various degrees of surface rough- 
ness (p), and two different subdominant order parameters. 



FIG. 4. Current density (normalized by Landau's critical 
current, j_t) at a (110) surface as function of the distance from 
the surface. Numerical results are presented for Tc^jTc ~ 0.3, 
T/Tc = 0.1, various degrees of surface roughness (p), and two 
different subdominant order parameters. Surface roughness 
is described by Ochchinnikov's model. 

The current resulting from a subdominant order pa- 
rameter of type s is an order of magnitude larger than 
from dxy, and much less sensitive to surface roughness. 
Important differences between surface states d.j.-2_y-2 + is 
and dx2-y2 -\-idxy are expected if the surface currents have 
to be matched along a closed surface consisting of, say, a 
(110), (1-10), (-1-10) and (-110) surface. The differences 
follow from the different symmetries of the two states. 
For instance, the state dr^-2_y-2 -\- idxy is invariant under 
the combined operation of time reversal and reflection 
X — >• —X (or y —>■ —y). As a consequence, the surface cur- 
rents of the state dx2-y2 +idxy can be matched smoothely 
in the four corners to give a closed current loop. On the 
other hand, the state dj.2^y2 ± is is invariant under the 
reflections a; ^ — x and y —>■ —y^ and the naive choice of 
a surface order parameter, namely the same state (e.g. 
dx2-y2 +is) on the (110), (1-10), (-1-10) and (-110) sur- 
face leads to surface currents with sinks (sources) at the 
corners. In order to have a single closed current loop for 
dg.2-y2 ± is-order, the surface states at adjacent surfaces 
have to be time-reversed. This requires at each corner a 
"defect" in the s-component at which dx2^y2 ±is changes 
into dx2-y2 =p is. This defect covers an area ~ ^q in the 
X — y plane, and thus costs little energy compared to 
the energy gained by forming the surface phases. The 
symmetry of the stable surface phase is determined for 
macroscopic surfaces by the channel with maximum sub- 
dominant Tc. Measurements of the surface transition 
temperature and the symmetry of the surface phase |25| ] 
provide new insights into the pairing interaction and thus 
into the mechanism of superconductivity in the materials 
of interest. 
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